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Abstract 

The infrared freezing of observables is known to hold at fixed orders of perturbative QCD if the 
Minkowskian quantities are defined through the analytic continuation from the Euclidean region. 
In a recent paper l] it is claimed that infrared freezing can be proved also for Borel resummed all- 
orders quantities in perturbative QCD. In the present paper we obtain the Minkowskian quantities 

CD : 

by the analytic continuation of the all-orders Euclidean amplitudes expressed in terms of the inverse 

n 

Mellin transform of the corresponding Borel functions |2j]. Our result shows that if the principle of 
analytic continuation is preserved in Borel-type resummations, the Minkowskian quantities exhibit 



a divergent increase in the infrared regime, which contradicts the claim made in yj. We discuss the 
arguments given in Q and show that the special redefinition of Borel summation at low energies 
adopted there does not reproduce the lowest order result obtained by analytic continuation. 
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I. INTRODUCTION 



Since the advent of QCD it was realized that the application of the renormalization- 
group improved perturbation theory is natural in the deep Euclidean region, where the 
running coupling is small and the physical hadronic thresholds are absent. The application 
of perturbative QCD for physical observables defined as Minkowskian quantities requires 
the analytic continuation from the spacelike to the timelike axis of the complex momentum 
plane. At high energies, the analytic continuation of the strong running coupling a(—s) from 
the Euclidean region s < to the Minkowskian region s > can be expanded in powers 
of l/ln(s/A 2 ). So, in the asymptotic region the expansion parameter is the same on the 
spacelike axis and the timelike one. At lower energies, however, one must take into account 
the finite terms appearing from the analytic continuation of lnf— s/A 2 ) — > ln(s/A 2 ) — m. 
The problem was investigated in the early 80' by several authors 3(, |^|, who tried to identify 
the most natural parameter for the perturbative QCD expansions of timelike observables. 
In j^] the authors compare the expansion parameters a(s),|a(— s)\ and Rea(— s) for s > 
and notice that \a{— s)\ seems suitable since it remains finite in the Landau region s < A 2 . 
However, the choice of the modulus |a(— s)\ as expansion parameter does not absorb all the 
7r 2 factors which arise from the analytic continuation, as shown by Radyushkin 4|, who 
derived explicit formulae for the timelike observables to every finite order of perturbative 
series. The analytic continuation was subsequently applied in the perturbative calculation 
of Minkowskian quantities |fj and in phenomenological analyses of inclusive observables like 
the rates of the processes e + e~ — > hadrons and r — > hadrons, using either low orders of 
perturbation theory [(|- [11 or resummations based on the Borel method jjj-jllj]. 

While for a long time the applications of perturbative QCD in the region < s < A 2 
were not considered reliable, the interest in the low energies increased when it was realized 
that some Minkowskian quantities, obtained in a consistent way by analytic continuation, 
remain finite in the timelike infrared limit s — ■> 0. This property called "infrared freezing", 
was shown to hold in every finite order of perturbation theory l^j], and is actually put 
on the basis of the so-called "analytic perturbation theory". In this approach 12 1, |l3j], the 
perturbative expansions of the Minkowskian observables are defined with a regular effective 
coupling, and the Euclidean quantities are obtained thereof by means of dispersion relations 
known to be valid in QCD under plausible assumptions. 
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One may ask whether the infrared freezing is only a feature of the finite order QCD 
expansions or it survives beyond finite orders. This is a nontrivial question, especially since 
the QCD perturbative series is known to be divergent. In using the Borel summation of 
the QCD perturbative series in the leading-/^ approximation, the authors conclude that the 
infrared finite limit of the Minkowskian observables is valid also at all orders in perturbative 
QCD. Since the perturbative series of QCD is ambiguous, it is not impossible, in principle, to 
implement a desired property by a suitable summation prescription. It is however natural to 
require that the procedure respects the principle applied to finite orders, which in the present 
case is the analytic continuation. In the arguments given in [J this principle is abandoned at 
some stage. The reason is that the authors use a Borel representation expressed as an infinite 
series of renormalons in the large-/?o approximation, which does not display the dependence 
on the momentum in a transparent way. So the question of what is the infrared limit of the 
Minkowskian quantities when defined in a consistent way by analytic continuation from the 
deep Euclidean region, as is done in the case of fixed orders, remains open. In the present 
paper we address this question. 

To this end, we choose an alternative representation of the Borel-summed Euclidean 
quantities, derived in Q], which is more convenient for the analytic continuation since it 
explicitly displays momentum dependence. A remarkable merit of this approach is that we 
do not need to represent the Minkowskian quantity in terms of any expansion parameter 
|a(s)|, a(\s\) or Rea(s) (as in IJ), assuming only that the quantity admits certain integral 
representations as discussed below in Section III. Note that the same technique was applied 
in [3] for the analytic continuation in the coupling plane, leading to results consistent with 
those obtained in 16j ]. 

As in jjj, we choose as Euclidean quantity the Adler function in massless QCD and as 
Minkowskian quantity the spectral function of the polarization function. In the next section 
we briefly review the analytic continuation from the Euclidean to the Minkowskian region 
of fixed-order perturbative expansions in QCD, stressing upon the fact that a consistent 
analytic continuation is free of ambiguities. In Section III we perform the similar analytic 
continuation of the whole Borel- resummed Adler function, written in a compact form in , 
which displays the energy dependence in an explicit way. In this Section we treat in detail 
the one loop coupling. The situation beyond one-loop is discussed briefly Section IV, where 
we show that our conclusion about the infrared behaviour of the Minkowskian observables 
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remains valid also in this case. We use the analytic expression of the two loop coupling 
derived recently in 17 1, working in the assumption, true in the real- world QCD, that 
the Euclidian coupling is not causal. In Section V we review the Borel summation presented 
in and show that it does not reproduce correctly the lowest order result obtained by 
analytic continuation in the infrared limit. 



II. ANALYTIC CONTINUATION OF FIXED-ORDER PERTURB ATIVE EX- 
PANSIONS 

We consider the Adler function in massless QCD defined as 

where n(s) is the correlation function of two vector currents. The function H(s) can be 
obtained from T>(s) by logarithmic integration: 

s 

n(s) — k — ln(-s) - J d ln(-s') V(s') , (2) 

where k is a constant and the integration is along a contour in the complex plane which 
starts at a fixed point and ends at s, without crossing the singularities of the integrand. 
This definition is consistent with asymptotic freedom and the general properties of the 
QCD Green functions. Causality and unitarity imply that n(s) and T>(s) are real analytic 
functions in the complex s plane (i.e. IT(s*) = n*(s) and T>(s*) = V*(s)), cut along the 
positive real axis from the threshold for hadron production at s = 4m^ to infinity. Along 
the cut, the Minkowskian quantity of interest is related to the spectral function ImiT(s + ie) 
by 

K{s) = -ImU(s + ie) - 1. (3) 
n 7r 
Following we consider the renormalization group improved truncated expansion of the 

Adler function in perturbative QCD 

V (N) {s) = a( _ s ) + J2d n a n+ \-s) , (4) 



n>l 



with the one-loop coupling defined as 



where A is the QCD scale parameter and (3q = (11N C — 2n/)/12 is the first coefficient of the 
(3 function (we follow in general the notations in Q], except for using (3q = b/2 instead of b). 
In our analysis we shall assume that /3q is positive, which means that infrared freezing does 
not hold for the Euclidian quantities like D(s) for s < 0. The first coefficients in (jljl. d n , 
n < 3, were calculated in jsj]. 

Using (J2J) we obtain the polarization amplitude as: 



U {N \s) = ife-ln(-s) 



In ln(-s/A 2 
Wo 



N 



n>l 



Jo, 



n+1 



nln n (-s/A 2 ) ' 



(6) 



We recall that the above expressions are derived in the deep Euclidean region s < —A 2 or, 
more generally, for complex values of s, with \s\ > A 2 . In this region the expressions are 
consistent with the general properties derived from field theory, which require that II(s) and 
T>(s) must be real for s < 0. The analytic continuations of (@J) and © to low values of \s\ 
contain however unphysical singularities on the space-like axis, which are absent from the 
exact amplitudes: T)( N \s) has a Landau pole at s — —A 2 , and n^(s) has a Landau cut 
along the interval —A 2 < s < 0. Here we are interested in the imaginary part of n^(s) on 
the upper edge of the timelike axis s > 0. Using (jBJ), the spectral function (J3J) is obtained, 
at finite orders, as 

N 

n( N \s)=A 1 (s) + Y,d n A n+1 (s), 



(7) 



where 



n>l 



Ms) 

AJs) 



1 



arctan(7r/5oa(s)) + 7r#(A 2 — s) 



Im 



(1 — inPoa(s)) 



l-n 



7T/?o Ti—l 

with a(s) defined in (j5)i. We note that in the first Eq. 



n > 1 



(8) 



arctan denotes the standard 



function defined in the interval (— tt/2, +7r/2), with arctan(O) = 0, and the term #(A 2 — s) 
accounts for the fact that the real part of ln(s/A 2 ) becomes negative when s < A 2 . Indeed, 
writing 

ln(-s/A 2 ) = ln(s/A 2 ) - in = yV^s/A 2 ) + tt 2 exp[i<j)} (9) 

for s positive and above the cut, one can see that the phase is continuous at s — A 2 , where 
it passes to the second quadrant (as shown in the first Eq. (JHJ) may be written also as 
A 1 (s) = l/TT(3 arccos[L/VL 2 + tt 2 ], where L = lns/A 2 = l/(/5 a(s))). 
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We mention that in some papers the Minkowskian quantity TZ(s) for s > A 2 is defined in 
terms of the Adler function through an integral along an open contour which ends at s ± ie. 
Usually, this contour is chosen as the circle of radius | s | centered at the origin, since in this 
case the integrals of the finite order expansions can be done analytically j^], While this 
procedure is suited for s much larger than A 2 , for points close to A 2 the result is sensitive 
to small deformations of the integration contour, which may or may not include the Landau 
pole. In y this ambiguity is solved by an ad-hoc choice of the branch of the arctan function 
which appears after integration, so as to lead to infrared freezing for finite order expansions. 
We stress that the procedure of calculating the discontinuity of the polarization function 
applied in the above Eqs. ©-© is free of such ambiguities. 

As was mentioned in the Introduction, in applications at large energies one expands the 
functions A n (s) in powers of the small coupling a(s) defined in (jSj). This gives for TZ(s) ^] 

K(s) ~ a(s) + d ia 2 (s) + (d 3 - di^p) a 3 (s) .... (10) 

The approximate expansion of TZ thus obtained can in no way be used at low energies, since 
a(s) becomes infinite at s = A 2 . On the other hand, as Eqs. (J7|) and imply, the TZ^ N \s) 
are regular for all s, including s = A 2 , and have a finite, universal infrared limit 



fcM(0) = i, (11) 
Po 



for N any positive integer. 



III. ANALYTIC CONTINUATION OF THE BOREL SUMMED AMPLITUDE 

The perturbation expansion (@J) of T>(s) in powers of the renormalized coupling a(—s) 
is known to be neither convergent nor Borel summable. We consider the Borel transform 
Bd{u) defined in the standard way in terms of the perturbative coefficients d n of V: 

Bo{u) = ± d ^(^) n , rf = l. (12) 

n=0 n - \P0j 

From the n\ large order growth of d n it is known that B^iuVhas singularities (ultraviolet 
and infrared renormalons) on the real axis of the w-plane For the Adler function, 

the ultraviolet renormalons are placed along the range u < U\, U\ = — 1 and the infrared 
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renormalons along u > 112, «2 = 2. Due to the infrared renormalons, the usual Borel-Laplace 
integral is not well-defined and requires an integration prescription. Defining 

oo±ie 

£>(±)( s ) = — / e -«/(A>«(-«)) b d {u) du = — lim / e -*/V>o*(-»)) B D {u) du , (13) 
Po J Po ^° J 

C± 0±ie 

one can adopt as prescription, for each value of a(— s) with Rea(— s) > 0, either T)( + \s) or 
or a linear combination of them, with coefficients £ and 1 — £ such as correctly to 
reproduce the known perturbative (asymptotic) expansion (J3J) of T>(s) (obtained by trun- 



1 these (and many other) 
the same asymptotic 



eating the Taylor expansion (J 12)) at a finite order N). Note that a 
integration prescriptions have, according to a theorem by Watson 
series in powers of a(—s) and therefore possess the same, original perturbative expansion. 

Once a prescription is adopted, one has a well-defined function, different prescriptions 
yielding different functions with different properties. In the present work we use, as in Q], 
the principal value (PV) prescription 

V(s)= 1 -[VM(s)+V^(s)}. (14) 

As shown in this prescription gives real values along the space-like axis outside the 
Landau region, which is consistent with the general analyticity requirements imposed by 
causality and unitarity. Moreover, we work in the ^/-scheme, where all the exponential 
dependence in the Laplace integrals (|T5|l is absorbed in the running coupling, and denote 
by Ay the corresponding QCD scale parameter. 

Our purpose is to obtain the Minkowskian qauantity 1Z by analytically continuing the 
Euclidean Borel-summed Adler function (114)1 . To this end it is convenient to use a repre- 
sentation of the Borel function Bp in terms of its inverse Mellin transform wb defined as 



21] 



2m 

uo—ioo 

The inverse relation 



uo+ioo 

wd{t) = — I duB D {u)r u - x . (15) 



B D (u) = jdrw D (r)r' u , (16) 

Q 

defines the function Bo{u) in a strip parallel to the imaginary axis with u\ < Keu < U2- 
The relations (JT5j) and (jT5j) are valid if the following L 2 condition holds 22]: 

uo+ico 

— j du\B D (u)\ 2 <oo, (17) 

Z7TZ J 



UQ—tOG 



where Un G [ui,u 2 }. The function w_d(t) was calculated in 

21] 

in the large-/3o approximation 



24j . where it has different analytic expressions, which we denote by and u> 



3 (» 



depending on whether r is less or greater than 1, respectively: 



W D {T) = Wp '{T), < T < 1 

w D {r) = w£\t), r > 1. (18) 

As discussed in jjj, one expects the inverse Mellin transform wp to have different expressions 
for r < 1 and r > 1 in general, also beyond the leading /^-approximation. Indeed, is 
given by a sum over the residua of the infrared renormalons, while is calculated in terms 
of the residua of the ultraviolet renormalons, and there are no reasons to expect these two 
contributions to be equal. In the large-/?o approximation, the expressions of the functions 



,«) 



2l|: 



w 



«), 

D I 



r ^--lnrJ + (1 + r) [L 2 (-t) + In rln(l + r) 
l + lnr+ Q + ^lnr^ + (l + r)[L 2 (-r- 1 )-lnrln(l + r- 1 ; 



(19) 



where L 2 (x) = — fn — lnfl — t) is the Euler dilogarithm. 

As noticed in |2l|, the function wd(t) defined in (fTHj) is continuous together with its first 
three derivatives, and satifies the normalization condition: 



J w D {r) dr = 1 . 



(20) 



On the other hand, Eqs. (|T9"j) define two independent functions, and which are 

analytic in the whole r complex plane except for logarithmic branch-points. These functions 
are not bounded everywhere: the function (r) is unbounded for r > 1 (growing at infinity 
like rln 2 r), while (r) grows like lnr/r for r — > 0. This behavior is seen in Figure 1, 
where we represent the function wd(t) defined by (fTHj). together with w^\t) for r > 1 
and w^\t) for r < 1. The same figure shows also the real parts of the functions w^\r) 
and w^\t) for r < 0, where they become complex. As we will see below, the growth of 
Hew^\r) for r — > — oo will have important consequences for the problem investigated in 
the present work. 

Following the technique described in detail in |2], we shall express the function T>(s) for 
complex values of s, with Rea(— s) > 0, in terms of the inverse Mellin transform €)£>. Using 




/ / 



w (t) for x >0 
Re w(<\x) for x <0 
wW(i) for x >1 
Re W^(t) for x <0 
w^(x) for 0<x <1 



-1 



FIG. 1: The function wd(t) denned by Q18|) in the large (3q approximation (i.e., with Wq(t) 
and wffi (t) given by (|19j0 is represented by a solid line. We separately display also the function 
ujp(r) for t > 1 (dash-dot-dotted line) and for r < (dashed line), and the function w^\t) for 
r < 1 (dash-dotted line and dotted line). For negative values of r the curves represent the real 
parts of the corresponding functions. The dashed line is unbounded for t tending to — oo, while 
the dash-dot-dotted line grows unboundedly for r increasing. The dotted and dash-dotted lines 
are unbounded for r approaching zero from either side. 

Eqs. (fT5|) and (fT^) as starting points, we rotate the integration contours C± in the complex 
w-plane up to a line parallel to the imaginary axis where the representation (J 16)) of Bp holds 
and can be inserted into the Borel integral. If the integrals are convergent, we can reverse 
the order of integration upon u and r, and perform first the integral upon the variable u, 
which can be done exactly. As explained in 3|, when s is in the upper half of the complex 
plane, the contour C + can be rotated towards the positive imaginary axis in the w-plane since 
the integrals remain convergent, while for the integral along the contour C_ it is necessary to 
first cross the real positive axis of the w-plane, picking up contribution of the residua of the 
corresponding singularities, i.e. the infrared renormalons. When s is in the lower half of the 
complex plane, convergence is achieved if the contours are be rotated towards the negative 
imaginary axis in the w-plane, and the roles of the contours C + and C_ are reversed. This 
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gives different expressions for V(s) in the upper/lower semiplanes of the s plane: 

V{s) = 1 /dr-^L + J ( ^ I ,^ '(-A? /,). im.s < 0. (21 ) 



/?o i ln(-r S /A y ) /3 





) 4 <f 











where the first terms are given by the integration with respect to u, and the last terms are 
produced by the residua of the infrared renormalons picked up by crossing of the positive axis. 
We recall that the expressions (|2Tjl were obtained by using the Principal Value prescription 

©• 

The corresponding expression for the polarization function H(s) can be obtained by in- 
serting the above result into the definition (J2J). This gives 



oo / \ . s 



n(s) =*- ln (i) -s/ dT ™ DWlnln (-i) ±f k l AH - s) (-tK<- a ^>- 

(22) 

with the contour in the last integral specified below Eq. (J2J) and the +/— sign corresponding 
to Ims > 0/Ims < 0, respectively The Borel-summed expression f!22j) . obtained for \s\ > 
Ay, can be analytically continued into the whole complex plane. We note that IT(s) is 
holomorphic for complex values of s and satisfies the reality condition II(s*) = n*(s). On the 
real axis, this function can have singularities manifested as discontinuities of the imaginary 
part. The unphysical singularities in the spacelike region —Ay < s < were discussed in 



detail in 



]. Here we are interested in the spectral function for s > 0. A straightforward 



calculation gives: 



A v /s 



1 / 7T \ 71 

Imll(s + ie) = 7i + — / drw D (r) arctan — -—5— + — / &tw d (t) 

Po{ \M r s/ A v)J Po J 



+ -^Re J driv^\r). (23) 



A. 

-Al/s 

We note that the last term in the first line was obtained by means of the relation 



Im 



lnln rS 



arctan I - - — - ) + n 9(A 2 V - rs) , (24) 
\\n(rs/A v ) 



already applied (for r = 1) in deriving the first relation (JHJ), and the last term in (J2*3*j) is 
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produced to the last term in Using (jSJ) we write also 1Z as 

oo . , A f-7 S 

n(s) = ±-Jdrw D (r) arctan (i^ry) + ^ / dr^r) + -Re / dr^(r). 

_A^/s 

(25) 

It is easy to check that this expression is continuous for all s > 0, including the point s = A v . 
We consider now the limit of this expression for s — > 0, i.e. A v /s — > oo. Since the first 
integral tends to zero (recall the comment below Eq. (jHJ) about arctan) and wd{t) satisfies 
the normalisation relation (J2*U|) . we obtain: 

1 °r 

H(0) = — + lim Re / dr^ <} (r). (26) 

fa J 

-A*, /a 

The first term coincides with the result ([lip, but we have now an additional term which 
involves the values of Wq(t) for negative r. Moreover, for small values of s, the integral 
involves arguments r — > — oo, where Rew^ (t) is unbounded (see Figure 1). Using ([19)1 
(which is valid in the large /3q limit) it is easy to check that the last integral in (|26J) diverges 
like ln 2 s/s 2 for s — > 0. This result disproves the statement made in |l|: the infrared limit 
of the Borel-summed Minkowskian observable lZ(s) obtained by analytic continuation from 
the Euclidian region does not reproduce the infrared freezing observed in the finite orders, 
and moreover displays an unphysical divergence. This behavior is crucially determined by 
the summation of an infinity of terms of the series. 

In order to understand the transition from the fixed orders to the resummed quantity, it 
is useful to apply the above formalism to the truncated perturbative expansion, when the 



Borel transform Bd(u) defined in (|T2*j) reduces to a polynomial B D 



(AO, 



u 



-/3o + 2!/3 2+ '-- + iV!/^ 



In this case the Laplace-Borel transform is well defined on the cuts. A straightforward 
calculation gives 

1 °r N 

/ e-/</W-)) B^\u) du=Y, d n a n+1 (-s) , d = 1, (28) 



n=0 



i.e. the finite order expansion (jlj). We want to check whether this expansion, along with 
the expansion (J7|) of TZ, are reproduced in the inverse Mellin formalism. 
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Clearly, when 5^ is a polynomial the condition (JT7j) is not satisfied, so we expect the 
function wd to be a generalized function (a distribution). In order to calculate it, we consider 
the alternative distribution function Wd(t), introduced in [icl ]: 

u +ioo 

W D {r) = ±-. du^r^. (29) 
2m J sm7ru 

UQ~iOQ 

When Bd(u) is a polynomial, the ratio Bd(u)/ sin7ru satisfies the condition fTTj) . which 
ensures the existence of the function Wd(t). For instance, for Bd(u) = 1 a straightforward 
calculation gives 

W D (T) = - n ^—r. (30) 

. . . . 7r(l + T) 

On the other hand, as shown in [lOj], [27[ (see also Eq. (A. 10) of j2j), the connection between 
the Mellin transforms Wd(t) and wd(t) is 

ir„(r) = iL^>, (31) 

7i J x + r 



where r can take arbitrary values, except for real negatives. By comparing (J30|) and (J31j) it 



follows that, for B K D '(u 



(0) 

K aj — x, 

(0) 



w%>(x) =8(l-x). (32) 



A straightforward calculation shows that at each finite order the function wd(x) is repre- 
sented in terms of the distribution 5(1 — x) and its derivatives. For instance, the inverse 
Mellin transform of B^ defined ()27|). for iV = 3, is 

^ 3) (x) = 5(l-x)-^5'(l-x)+^r(l-x)+y(l-x)]-^[r(l-x)+35''(l-x)+5'(l-x)] . 

P0 2\p 6\/J 

(33) 

Such a representation is not unique: except for the first two terms which remain the same, 
the higher terms can bewritten equivalently as the product of the n-th derivative of 6(1 — x) 
with a polynomial of degree n — 1. For instance, w D (x) in ()33|) can be expressed in the form 

w%\x) =6(l-x)-^5\l-x) + 7 ^- 2 [(x + l)/2 8"(l-x)]-^[(x 2 +Ax+l)/6 8"\l-x)\. 

(34) 

It is easy to check that the different expressions ()33|) and give the same result for the 
quantities of interest T> and 1Z. An immediate consequence of these expressions is that, for 
finite orders, Wd < ( t ) = 0- Inserting tD® into the relation ()21|) . which expresses the Adler 
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function in terms of the inverse Mellin transform, we obtain by a straightforward calculation 



1 T r°° / \ r^ 2 / s 

K (0) (s) = — — / 5(1 — x) arctan ( 7r/(ln x + ln(s/A 2 )) ) dx + 7r / 6(1 - x)dx 
TVPn [Jo v ' JO 

arctan(7r/? a(s)) + ty9(A 2 — s) , (36) 



(3 \n(-s/A 2 v ) \(3 ln(-s/A 2 v )J \/3 Hs/K)) \(3 ln{-s/A 2 v ) J ' 

(35) 

which coincides with the first terms in the expansion (J3J). 

Let us insert also into the resummed expression (|25[) of 7?.. It is easy to see that the 
first term Wp\r) = 6(1 — r) contributes both to the first and the second integrals in 
giving the result 

/ ^( 1 — .rl arctan ( tt I ("In x 4- In f.s / \ z ) ) ) dx A- it I 

1 

7T/3 - 

which coincides with the function Ai(s) defined in (JHJ) and satisfies the property of infrared 
freezing. The higher terms in (}3*3*j) contribute only to the first integral in (|25|h reproducing 
the terms in expression (|25|) . For instance, inserting the second term of (|33|) in (|25|) one has 

J 6'(1 — x) arctan [re/ (Inx + ln(s/A 2 ))J dx 

with A 2 (s) defined in (|SJ). So the formalism of inverse Mellin transform reproduces the finite 
order expansions which are consistent with the property of infrared freezing But the 

summation of the whole series leads to a different result. The discussion in this Section 
reveals the difference between the finite orders and the summed expression: it resides in the 
function w^(t), which is zero at each finite order but is nonvanishing when an infinity of 
terms are summed and the infrared renormalons show up. 



n(3 2 



dl l + a 2 PW =dlMs) > (3?) 



IV. BEYOND THE ONE-LOOP COUPLING 

Up to now we restricted the discussion to the one- loop coupling (J3J. We show now that 
the same conclusion is valid beyond this approximation. If the one loop coupling is not 
inserted into p3jl . it is easy to see that the two integrals in (|2*2*J) write in general 

s oo s 

n(s)~--=- /dln(-s) fdr- — ^4 ±^ /dln(-s)e~5^F^« ( e ~s^) . (38) 

fo J jL lnr + M^) P° J V J 

As shown recently jl7[ , [18[ the solution of the two-loop /3-function equation can be written 

analytically in closed form as 

1 1 / S \-Po/c 

a{s) = - c[l + W( Z (s))] ^ zis) = -e{T*) ' (39) 
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where c = (153 — 19n/)/24/3 is the second universal beta-function coefficient and W(z) is 
the Lambert function defined implicitly by W(z)e 



W{z) 



We work in the condition 



c > 0, valid in real- world QCD. Then, as shown in [l^j], the physical branch of the Lambert 
function in (J39|) is and the coupling does not freeze in the spacelike region. 
Denoting W(z) = W-\(z(— s)) we obtain from (1391) : 

s \ -/V c 



A 2 



-W(z)e 



W(z)+1 



and 



ln(^)=-^(ln(-l) + ln^) + l) (40) 

c c c 1 +W 

dln(-s) = -— dlnz(-s) = -— (dlnW + dW) = -- — - — dW. (41) 

We now insert the two loop coupling (}3*9~j) into (}3*Hj) . and notice that the first integral can be 
performed by making the change of variable (|41j). Using the relation W'(z) = W(z)/z(l + 
W(z)) (to be obtained from We w = z by differentiating the logarithm) we write after a 
straightforward calculation the first integral in ()38|) as 



1 °° 

11 / dr «^( r ) 

^° o 



• (42) 



c / 1-^lnr l + W(z)-^lnr 
In order to evaluate the limit s — ► 0, we use the asymptotic expansion j3] 

W(z) ~lnz-lnlnz (43) 

valid for both z — > and z — > oo. Then we obtain from the first logarithm in (J42j) the limit 

oo 

Mmn(a) ~ -- /dr^ D (r) Lain (^) , (44) 
^° o 

while the second term in (f4*2"|) vanishes. The expression (jH| coincides with the first integral 
in (j22J), whose imaginary part, as discussed in Eqs. (|2*3~j) - (j2fi)) . has the infrared limit l/3n. 
This result is consistent with the finite order calculations with the two loop coupling in |l|, 



26|. 



We turn now to the second integral in (|38|). It is easy to see that the leading term in 
the asymptotic behavior ()43)1 gives, for s — > 0, an expression identical to the second term 
in (j2SJ- The next to leading term in (|43jl introduces logarithmic corrections which do not 
change the singular behavior of the integrand for s — > 0, i.e. t — >• — oo. Thus, although the 
two loop coupling leads to a different behavior of the spectral function at nonvanishing s, 
the singularity in the infrared limit is dominated by the one loop coupling. 
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V. COMMENTS ON THE PROOF OF INFRARED FREEZING IN [1|| 

In Ref. jl| the authors discuss the infrared freezing beyond fixed-order expansions by 
using a representation of the Minkowskian quantity 1Z in terms of the Borel transform Bp 
of the Adler function. This representation is derived by inserting the Borel representation 
of the Adler function into the expression (JHJ) of the polarization function, which is then 
used in (J3J) to compute the quantity 71. Adopting the Principal Value prescription (JHJ), a 
straightforward calculation gives 

Q.Q: 

1 °° 

K(s) = -^PV / e -«/^» B D (u)du, (45) 



with the one-loop coupling a(s) defined in 

The integrals in (??) converge only for s > A 2 , when a(s) > 0. For small s, of interest 

for the infrared behavior, a(s) is negative and the standard Laplace-Borel integral, along 

n 

the positive real axis of the w-plane, diverges. In |l( the authors notice that for a(s) < 
a convergent integral is obtained by choosing as integration line the negative axis in the u- 
plane, instead of the positive one. Therefore, they define the Borel-summed 7Z{s) for s < A 2 
as: 

n(s) = -^PV / e- w /^°°( s » B D (u)du. (46) 

The PV prescription now regulates the ultraviolet renormalons along the negative real axis, 
while in the standard definition (??) the prescription regularizes the infrared renormalons. 
The authors of l| claim that by this redefinition of the Minkowskian quantity for s < A 2 , 
one recovers infrared freezing (fTTj) also beyond fixed-order perturbation expansions. 

In principle, a redefinition of the Borel integral as in (??) is not illegitimate: since we deal 
with functions which are not Borel summable, the choice of the prescription for the Borel 
integral is to a large extent arbitrary; as mentioned in section III., there are many different 
functions that have the same perturbative, divergent asymptotic expansion in powers of the 
coupling constant. Such a redefinition should however be motivated physically, not just by 
the fact that the original definition has lost mathematical sense in a region. 

Furthermore, when introducing a redefinition of the Borel integral, we must ensure that 
the fixed order expansions are reproduced when the Borel transform is expanded in a Taylor 
series. In [lj the authors claim this requirement to be fulfilled if one expands in powers of 
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u not the whole integrand in (??), but only the Borel transform Bd(u). However, it is easy 
to see that the prescription (??) fails to reproduce correctly the infrared freezing (jllj) for 
finite order expansions. 

Indeed, let us insert in (??) the truncated expansion (J27j) of B^(u) in powers of u. Then 
the Borel integral is regular, no prescription is required and we obtain from (??), for s < A 2 : 

— oo oo 

n W( 8 ) = 4- I e~ u ^<^ S -^B^(u)du = ~ / e -"/(/?o(-^)) s J^ B ^\-u)du, 

3q J TTU On J TTU 



(47) 

where a change of variable was performed in the last step. In the last integral the quantity 
—a(s) is positive for s < A 2 (recall the definition (JSJ) of the coupling), and we can easily 
perform the integration for each term of B^\—u) from (|27j) . A straightforward calculation 
shows that, for all the terms except the first one, the infrared limit of the integral (??) is 
zero, in agreement with the behavior of the functions A n (s) of (JHJ, with n > 1. For the first 
term Bd(u) = 1 (which was responsible for the infrared freezing of the truncated expansion 
in Section II) we apply the identity (41) of [l| which gives 

oo 

K®(s) = --=- / e-^-MH du = -4- arctan[vr/? (-a( S )] . (48) 

PO J TT/3 

In the infrared limit s — > 0, Eq. (jSJ) implies— a(s) — > through positive values, and from 

(??) we obtain TZ^(O) = 0, which is not consistent with the infrared limit of the function 

fl 

A 1 (s) defined in (JHJ), and with the relation (Jllj) . Therefore, the prescription adopted in Jl| 
for the Minkowskian quantity at s < A 2 fails to reproduce correctly the infrared freezing 
(I11J) of the truncated expansion. 



VI. CONCLUSIONS 

In a recent paper Q it is claimed that by using analytic continuation in the energy 
plane it is possible toprove the infrared freezing of Minkowskian quantities beyondfinite- 
orders in perturbative QCD. In the present work, we applied the technique of the inverse 
Mellin transform of the Borel function, developed in j2j], which givescompact expressions 
of the QCD amplitudes in the complex plane. By using the analytic continuation of these 
expressions into the Landau region, we calculated explicitly the spectral functions, as in the 
fixed order expansion. As in [lj we adopted the Principal Value prescription, and considered 
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as Euclidean quantity the Adler function in massless QCD. Our result, expressed in Eq. 
(|26|). contradicts the conclusion reached in the summation of higher orders in QCD 
leads to a divergent increase of the Minkowskian quantities in the infrared limit, if these 
are calculated by analytic continuation from the Euclidean region. The divergent infrared 
behavior arises explicitly from the summation of the infinite terms and is related to the 
infrared renormalons. Of course, one expects that in full QCD this divergent behavior will 
be compensated by a similar growth of remaining terms in the OPE, calculated with the 
same prescription. 

The difference between our results and those in [J is explained by the fact that the 
authors of [j| do not apply consistently the principle of analytic continuation, applied at 
finite orders. Lacking a compact expression of the Adler function, as the one provided by the 
inverse Mellin transform used by us, the authors make the analytic continuation of the Borel 
representation itself, which is valid only outside the Landau region. Therefore, the Borel 
representation (??) of the Minkowskian quantity considered in |l[ converges only for s > A 2 , 
and is useless in the infrared limit. To reach this point the authors change the definition 
of the Borel integral. However, as we showed in Section V, this new prescription for the 
Borel summation of Minkowskian quantities below the Landau point fails to reproduce the 
infrared freezing of the truncated expansion . 
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